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For adiabatically and periodically manipulated dissipative quantum systems we derive, using Flo- 
quet theory, a simple Markovian master equation. Contrary to some previous works we explicitly 
take into account the time dependence of the Hamiltonian and, therefore, obtain a master equation 
with a time-dependent dissipative part. We illustrate our theory with two examples and compare 
our results with the previously proposed master equations. In particular, we consider the problem of 
Cooper pair pumping and demonstrate the inadequacy of the secular (rotating wave) approximation 
when calculating the pumped charge. The secular approximation producing a master equation of the 
Lindblad type approximates well the quantum state (density matrix) of the system, while to deter- 
mine the pumped charge a non-Lindblad master equation beyond the rotating wave approximation 
is necessary. 
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I. INTRODUCTION 

The hope to develop a quantum information process- 
ing machine has made a strong impact in several distinct 
branches of science, e.g., in theoretical and experimen- 
tal quantum physics, in computer science, and in elec- 
trical and chemical engineering. To succeed, it is essen- 
tial to achieve high-precision quantum transformations, 
which so far are hindered by two problems: decoherence 
(see Jj for notation) and gate control. While the former 
can be solved to some extend by proposing new types of 
less fragile qubits, or implementing decoherence free sub- 
spaces [21 [3] , the solution to the latter could be the use 
of adiabatic transformations @]. Such transformations 
have the advantage that they can be surprisingly robust 
to various errors of experimental parameters. 

This robustness of adiabatic gates, however, comes at 
the price of a significantly longer gate time, therefore 
giving the environment more time to destroy coherence. 
For this reason, it is important to understand the inter- 
play of adiabaticity and decoherence, which was studied 
in several publications in connection to Berry's phase [SI- 
IS]. In these, the authors usually postulated a Markovian 
master equation of Lindblad form. 

However, microscopic derivations of master equations 
show that its structure crucially depends on the Hamilto- 
nian of the system. Therefore, the time dependence of the 
Hamiltonian should be explicitly taken into account when 
deriving a master equation. This approach was already 
taken, e.g., in Refs. [TUHH] which studied the pumped 
charge in a Cooper pair sluice, as well as in Ref. [13] and 
Ref. [14 which studied hysteresis in molecular magnets 
and entanglement of two spin molecules, respectively. In 
this paper we are interested both in general properties of 
decoherence in adiabatic systems as well as in the con- 
crete problem of Cooper pair pumping. We also discuss 
the applicability of the secular approximation and show 
that it leads to an incorrect pumped charge despite the 



fact that this approximation is very well justified if one is 
interested in the quantum state (density matrix) of the 
system. 

We wish to point out that there are two distinct 
notions of adiabaticity. First, in the quantum litera- 
ture one speaks of adiabatic behavior if a closed sys- 
tem approximately follows the instantaneous eigenstates 
of the Hamiltonian H{t), i.e. \ip{t)) « e*'^(*) \n{t)) with 
H{t) \n{t)) — En{t) \n{t)). Second, in the thermodynam- 
ical literature adiabaticity means that the system fol- 
lows the instantaneous thermal equilibrium state p{t) cx 
exp[~H{t) j (/csT)], where T is the temperature of the en- 
vironment. While both situations require a slow change 
of the Hamiltonian, in the quantum case the time scale is 
set by the eigenfrequencies of the Hamiltonian, whereas 
in the thermodynamic case the energy relaxation time 
(often referred to as T\ in the context of qubits) is the 
appropriate time scale. A master equation describing an 
adiabatically steered quantum system which is coupled 
to a thermal bath should therefore recover the thermo- 
dynamical adiabatic behavior in the respective limit. 

If Tr is much longer than the period T of a cycli- 
cally time-dependent Hamiltonian, then one would ex- 
pect that in the quantum adiabatic limit, the system ap- 
proximately follows a mixture of the instantaneous eigen- 
states p{t) fa J2m ^m{t) \rn{t)) {m{t)\ (or more precisely of 
Floquet states [E])- After a transient period the coeffi- 
cients Cn{t) approach constant values c„(cx)) which de- 
pend on some sort of time-averaged level splittings of 
the Hamiltonian as well as on the temperature T of the 
environment. This is the limit which was described in 
Ref. [ig. 

However, in the other limit, <C T, one expects 
that after some transient time the system follows the in- 
stantaneous thermal state. The coefficients then satisfy 
c,n{t) = c„(i) exp[{En{t)-Emit)/{kBT)] and hence cycli- 
cally depend on time. Our work represents an extension 
of Ref. [12] as both limits as well as the crossover between 



them can be described by our master equations Eq. ( 14 ) 
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and Eq. (17) 



In section [II] we introduce the Floquet formalism for 
cyclic Hamiltonians. We then modify the formalism 
to separate the fast dynamics which arise from the en- 
ergy splitting of the Hamiltonian from the slow dynamics 
which are due to the adiabatic time dependence of the 
Hamiltonian. This procedure directly points to an ap- 
proximation suitable for adiabatically driven, Markovian 
open quantum systems, which results in a Lindblad type 
master equation, but with time-dependent Lindblad op- 
erators. We find that decoherence is best described in 
the basis of the time-dependent Floquet states in the 
sense that in this basis one can differentiate between de- 
phasing (the environmental influence on the off-diagonal 
elements of the density operator) and relaxation (diago- 
nal elements). The decoherence rates generally depend 
on time to guarantee that in the thermodynamical adia- 
batic limit the system follows the instantaneous thermal 
state. It follows that a system coupled to a zero temper- 
ature environment does not relax into the instantaneous 
ground state of the Hamiltonian, but into the slightly dif- 
ferent Floquet ground state. We apply our theory in two 
simple examples. In one of these examples, the problem 
of Cooper pair pumping, the measurable observable is the 
transferred (pumped) charge. We show that although the 
Lindblad master equation leads to a very good approxi- 
mation for the density matrix of the system, it predicts 
a completely wrong transferred charge. To rectify the 
problem we refrain from the secular approximation and 
obtain a suitable master equation which is not of Lind- 
blad form. 



II. FLOQUET THEORY 

In this section, we first present some well-known re- 
sults of the Floquet theory [151 EU ■ We then modify this 
theory to allow for a separation of fast oscillations due to 
the energy of the states from the slow oscillations due to 
the periodic nature of the Hamiltonian. 

For a time-periodic Hamiltonian H{t) = H{t + T) 
there exist solutions of the Schrodinger equation of the 
form 



l*a(0) 



l*a(0): 



(1) 



where the Floquet modes = \'^a{i + T)) have the 

same period as the Hamiltonian, and are the quasi en- 
ergies. The quasi energies are only uniquely defined up to 
multiples of Ml with f2 = 2i:/T, because e""^* \'^a{t)) 
for all n G Z correspond to the same Floquet mode but 
with a shifted energy. 

Adiabatic Floquet theory. For our purposes, it is 
advantageous to separate the fast dynamics associated 
with the instantaneous eigenfrequencies of the Hamilto- 
nian from the slow dynamics associated with the time de- 
pendence of H{t). Such a separation is not provided by 
Eq. ([T]) because the phase factors in Eq. ([T]) oscillate with 



constant frequencies Eq. Therefore, if the eigenfrequen- 
cies depend on time, as is generally the case, the Floquet 
modes [^^(t)) also oscillate rapidly. A more convenient 
representation is provided by the states |(/)Q,(i)), which 
differ from the Floquet modes only by a time dependent 
phase factor: 

|*a(i)) - e-*-^"'^*'-^"(*')e-"^°*/'^|(/.c«(t)). (2) 

Here, Ea{t) — {ipa{t)\ H{t) {(pait)) is the instantaneous 
energy of the Floquet mode, and if a is chosen such 
that the states \4'a{t)} = \4'a{t + T)) are periodic in 
time. In fact, it can be easily shown that ipa — 
i dt {(l)a{t)\ \(j)ait)) is the non-adiabatic geometric 
phase (Aharonov-Anandan phase) [T7] associated with 
the state \(j)a{t)). 

Note that the ambiguity in the quasi energies of the 
original Floquet states is non-existent in the modified 
Floquet states. Therefore, we can now speak of a well 
defined Floquet ground state and Floquet excited states. 
Although the concept of a Floquet ground state is not 
very meaningful if the driving is fast, for adiabatic driving 
it will turn out to have the same importance as the energy 
ground state for non-driven systems. Indeed, the Floquet 
ground state approaches the instantaneous ground state 
in limit of T ^ oo. 

In the adiabatic limit, a system prepared in an en- 
ergy eigenstate will stay in the respective instantaneous 
eigenstate. Because the eigenstates are cyclic, they are 
identical to the modified Floquet states in a zeroth order 
adiabatic approximation. Therefore, if the time evolution 
is adiabatic, the following approximations are valid 



\Mt)) = M)) + o{A), 



(3) 



Here, \na{t)) and Ea(t) are the instantaneous eigen- 
states and energies of the Hamiltonian, respectively, (fa 
is the adiabatic geometric phase (Berry phase), and 

A — maxt^Q^^ jg ^Yie adiabatic parame- 

ter. These modified Floquet modes have the useful prop- 
erty to vary only slowly in time, such that the expansion 



oo 



\Ca,k) 



(4) 



k— — oo 



can be restricted to \k\n ^ min^^^.t \Ea{t) — Eis{t)\. This 
property will turn out to be useful in the derivation of a 
master equation. 

Super adiabatic expansion. Later on we will find 
that the system decoheres into the Floquet states. That 
is why we want to have a closer look at the relation of the 
Floquet states to the instantaneous energy eigenstates. 
The statement that an adiabatically evolving system fol- 
lows the instantaneous eigenstates of the Hamiltonian 
H{t) is only correct to zeroth order of the adiabatic pa- 
rameter. It is well known 1181 that if the Hamiltonian 
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is an analytic function in time, then the system follows 
more closely the eigenstates of the operator 



H ^ H 



(5) 



where we omitted the time dependence of the oper- 
ators for shorter notation and used a dot to indi- 
cate the time derivative. The operator Uo{t) is de- 
fined by \n{t)) ^ U^{t)\n{0)). Similarly, Ui{t) is de- 
fined by \ni{t)) — ul{t)\ni{0)), where \ni{t)) are the 
eigenvectors of the transformed Hamiltonian Hi{t) = 
Uo{t)H{t)ul{t) + iUf){t)ul{t). The first term in Eq. ^ 
is of zeroth order in the adiabatic parameter, the second 
of first order, and so on. 

One should note that for any given adiabatic parame- 
ter, the series Eq. (5]) does not (except for some special 
cases) converge [TH . The smaller the adiabatic param- 
eter, the longer the series improves before it starts to 
diverge. If the series does converge, then the system will 
exactly follow the eigenstates of H which are then the ex- 
act Floquet states. Even if the series does not converge, 
in an approximation which improves exponentially with 
the adiabatic parameter, one can use the first few sum- 
mands of Eq. ([5| to obtain the Floquet states. Thus, by 
using the Floquet modes our theory implicitely takes into 
account terms of higher order in the adiabatic parameter. 



III. MASTER EQUATION 

In this section, we derive a master equation in the mod- 
ified Floquet basis. Some properties of this basis will 
lead to approximations which are appropriate for adia- 
batically evolving systems. 

We assume that the total system evolves under the 
influence of the Hamiltonian 



HAB{t) = HA{t) + HB+A®B, 



(6) 



where Ha and Hb are the Hamiltonians of the system 
of interest and the environment, respectively, and A® B 
is their coupling with A acting on the Hilbert space of 
the system of interest and B acting on the environments 
Hilbert space. We assume an initially uncorrelated state 
p(0) = pa(0) ® Pb where pb describes the thermal equi- 
librium of the environment. Without loss of generality 
we assume Tt{pbB) = 0, and using the Born-Markov 
approximation we arrive at the well known interaction 
picture master equation [19] 



/•OO 

p^Ait) = / dTTTB[B{T)BpB] 

Jo 



(7) 



X [A{t-T)pA{t)A{t) - A{t)A{t-T)pA{t)] + c.c, 

where A{t) and B{t) denote operators in the interaction 
picture. We note that the unitary evolution operator 
used to transform A into the interaction picture may be 
written conveniently as 



C/^(t) = ^|vI'„(<))(0„| 



(8) 



where we introduce the notation \(j)a) = \(j)a{0))- We use 
this form of the evolution operator to write 

A{t) = u{{t)AUA{t) 

= 5^(^'„(t)|A|vI/„,(i))|0„)(0„,| 



aoi.'k 



(9) 



where we used the Fourier expansion 



Aaa\k — 



r 



6„(i)| A (0), (10) 



as well as 

UJaa'it) = ^«'(i)-^a(i) + (<^a' -<^a)/r (11) 

For most purposes, the angular frequencies uJaa' = {Ea' — 
Ea)[^+0{A)\ may be approximated by the instantaneous 
transition frequencies. 

We now substitute Eq. ^ into Eq. ^ to find the 
Bloch-Redfield type equation 

'{t) + k^]Aaa' ^kAp' p^l 

X e'-l^o d-*' ["3/3' (*')-"»„' (t')-(fc+/)n] 

X [\(l)a){(t)a'\PA{t) \(t>P'){(l)0\ 
- l0/3')('^/3|0a)(0a'|PA(i)] 

+ C.C., (12) 

where we define the one-sided Fourier (Laplace) trans- 
form of the environment correlation function 



r[uJaa'{t) + kH] — / dr e 



,i//_dt'K„,(t')+fcf2] 



TrB[BiT)BpB] 



/ dr e*['^-' (*)+'=^f^l^TrB[B(T)SpB] . (13) 

^0 



The latter form, which is more familiar from the theory of 
time independent systems, is valid if the eigenenergies do 
not change much on the time scale set by the correlation 
time of the environment. 

Using the standard Floquet modes, the master equa- 
tion ( 12 ) was also derived in the review of Grifoni and 



HanggFflS] ■ and served as the starting equation in |12| . 
where the authors proceeded by performing two rotating 
wave approximations (RWA). The first RWA is to ne- 
glect terms with k ^ I, and the second one is to keep 
only terms with uJaa' = The first RWA is valid, 

if g-»(fc+0f2t can be considered as fast rotating. In par- 
ticular for adiabatically evolving systems this might not 
always be the case, and therefore, we do not use this 
RWA in the following. 

We perform an approximation which is particularly 
suitable for adiabatic driving, namely T{uj + kfl) w r(a;). 
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The requirement of a sufficiently smooth correlation func- 
tion T{uj) is satisfied if the modified Floquet modes 
change slowly compared to the environment correlation 
time. 

To arrive at a master equation of Lindblad form and 
to guarantee complete positivity, we have to perform the 
second above described RWA, usually referred to as the 
secular approximation (although, as we will see later, 
there are some instances where the secular approximation 
leads to erroneous results and one has to forego complete 
positivity at all times). Noting that in general a RWA is 
valid if the exponential in the neglected term is rotating 
fast, we recognize from Eq. ( 12 1 that this RWA requires 
that (k — l)fl is not of the same order as ujaa' — i-^pp'- 
This can only be guaranteed if the Hamiltonian is varied 
adiabatically, and if the expansion in Eq. Q can indeed 
be restricted to only small values of k. Therefore, the use 
of the modified Floquet modes introduced in this paper is 
essential for this RWA to be valid, even in the adiabatic 
regime assumed here. 

Noting that Aaa',k vanishes unless kfl is much smaller 
than the transition frequency WQC['(i), we proceed by ne- 
glecting all terms of Eq. ( 12 ) with w^a' (t) ^ ojpp' (t). Us- 



-ikVtt 



A. 



('/'op) I A \(t)p(t)) and transforming 



back into the Schrodinger picture, Eq. \12\ becomes 

PA{t) = -i[H{t) + Hi^s{t),PA{t)] 

+ 7(0) Lo{t)pA{t)Ll{t) - ^ [Ll{t)Lo{t),pA{t)} 

+ 51 '^['^"'3(0] 

ap{t)pAmipit) - i {Lipit)La„it),pAit)] 

(14) 



with the Lindblad operators 

io(t) = ^(0a(i)|A|0„(t))|0„(i))(0„(i)|, (15) 

a 

Lo.p{t) = (0a(t)|A|0^(t))|0„(t))(0^(O|. (16) 

While the operator causes pure dephasing (i.e. it 
does not change the populations of the Floquet modes), 
the operators Lap are responsible for energy relax- 
ation. We further separated the imaginary and real 
part of the Laplace transformed environment correla- 
tion function r(a;) = ^7(w) + iR{uj) and found the 
Lamb shift contribution to the Hamiltonian i?Ls(i) = 
J:^^ R[u;ap m I (0a it)\A\Mm'\ ^ (O) W l to be di- 
agonal in the Floquet basis. 

Thus, taking into account the explicit time depen- 
dence of the Hamiltonian, we showed that the system 
does not decohere into the instantaneous eigenstates, 
but into the Floquet states, a result previously obtained 
by [12]. The two bases are equal in the adiabatic limit, 
but they differ for finite periods T as according to 



Eq. ([5]) the Floquet states implicitly include higher or- 
der adiabatic contributions. While pjj uses rates which 
are averaged over one period, we take into account the 
rates '-f[uJap{t)]\{(t)a{t)\ A \(j)p{t))\'^ at each instant of time, 
which might be important in particular for adiabatic 
driving. On the other hand, the time dependent rates 
were already obtained in [T31 [2] , where the authors as- 
sumed an infinite period 7" — > 00 and hence found that 
the system relaxes into a mixture of instantaneous energy 
eigenstates. Therefore, our theory generalizes previous 
works and reduces to the known results in the respective 
limits. 

Because of the time dependent rates of Eq. ( [l4| , we 
can not expect to find a steady state in the Floquet basis. 
Nevertheless, there is a quasi stationary state pA{t+T) = 
PA{t) with vanishing off-diagonals in the Floquet basis. 
If the temperature of the environment is zero, i.e. T{lj < 
0) = 0, then the system evolves into the Floquet ground 
state. 

In the application of our theory to Cooper pair pump- 
ing, we are advised by the well known fact [TUl HH HO] 
that the use of the secular approximation may lead to an 
error in the transferred charge, despite the fact that it 
leads to a very good approximation of the actual density 
operator (as explained in subsection IV B I. Therefore, we 
provide here the Schrodinger picture master equation ob- 
tained without the secular approximation. This master 
equation follows from Eq. (12 1 with r(cj + kil) w r(a;). 
Thus we obtain 



PA{t) 



with 



~t[HA{t),pA{t)] + K{t)pA{t)A^ +ApA{t)K^{t) 
-A^K{t)pA{t)- PA{t)K\t)A (17) 



K{t) = ^r[w„„,(0] (0c(i)| A \(^a'{t)) \(l>a{t)){(l>a'{t)\ . 



(18) 



This master equation includes the cross-coupling between 
diagonals and off-diagonals of the density operator. In 
what follows we will use both, Eq. (14) and Eq. (17 1, 
depending on the situation. 



IV. ILLUSTRATIVE EXAMPLES 

In this section we illustrate our theory with some sim- 
ple examples of two-level systems, which could poten- 
tially be tested in experiments. In the first, the Floquet 
states can be found analytically. It serves to demon- 
strate the time dependence of the energy relaxation and 
dephasing rates, as well as the relaxation into the Floquet 
ground state at zero temperature. The second example 
is slightly more complex as it has a time dependent level- 
splitting. Therefore, it exhibits a quasi stationary state 
where the populations in the Floquet basis depend on 
time. 
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A. Two level system with constant level— splitting 



We consider the Hamiltonian 



H{t) = \B\\a, 



1 



B±\cos{Vlt)ax + sin(rit)(Tj, 



(19) 



where the matrix form is written in the basis |t) , |i) 
of eigenvectors of a^- The adiabatic parameter is easily 
obtained to be = VlB±/{B\ + -B[|). This Hamiltonian 
describes, for example, a spin in a magnetic field which 
rotates around the z-axis. Because the eigenvalues of 
the Hamiltonian do not depend on time, the modified 
Floquet states equal the traditional Floquet states: 



|0+(t))=a|t>4 
\ci^At)) = -h\^) 

with real a and h satisfying 
a/h ^ 



iVtt 



he 

\- ae 



ID, 



25 I 



Bii -n + J{B\\ -n) 



Bl 



(20) 



(21) 



and + 6^ ~ 1. The corresponding Floquet energies are 
given by 



n±J{B\\ - ny 



Bl 



(22) 



and are equal to E± + ip±/T. 

We assume a coupling of the system to the environ- 
ment via Gx , which could be a fluctuating external mag- 
netic field in x-direction. Therefore, we find the Lindblad 



operators in the Schrodinger picture from Eq. ( 16 1 with 
A = cJx 

Lo = 2a6cos(m)(|</.+ (t))(0+(i)l"l'/'-(i))(0-(<)l), 



-int 



5^e^^" |0_(i))(0+(i)| ,(23) 



We see that there are times when dephasing is the domi- 
nant decoherence effect {Lq), and times when relaxation 
(L |_) is dominant. 

In Fig. [l| we show results of the numerical integra- 
tion of the master equation for an initial state \4>+)- 
In panel (a), the coupling to the environment is strong 
enough to force a thermal mixture of the Floquet 
ground and excited states within one period T, while 
in panel (b), several periods are needed for the system to 
evolve into the thermal mixture. The populations of the 
Floquet states approach a constant in the long time limit 
because the level splitting of the Hamiltonian does not 
depend on time. The time dependence of the relaxation 
rate is noticeable in both panels. 

At zero temperature, the system relaxes into the 
Floquet-ground state |(/)_(i)), rather than into the in- 
stantaneous ground state of the Hamiltonian. Note that 
if the time dependence of the Hamiltonian were not ex- 
plicitly taken into account in the derivation of the master 




FIG. 1: (Color online) Panel (a) shows the population of 
the Floquet excited state \(f)+{t)) during one adiabatic cycle. 
The damping 7(1:1; = (t+ — t-)/K\ = 0.7 Q,f? is strong com- 
pared to the period T, such that the quasi stationary state 
is closely approached within one period. In panel (b) the 
damping ^[uj = (e+ — = 0.07 is weak such that 

the quasi stationary state is only reached after many periods. 
While black lines shows the populations of the Floquet ex- 
cited state |0+(t)), blue lines are the populations of the state 
It). The gray dashed line is the population of |t) if the sys- 
tem would be in the instantaneous ground state. Parameters 
are chosen as B\\/h = 2 0,, B±/h = 5 f2, and temperature 
kBT/{2nh) = 00 (solid), 5 fl (dashed), 10 Q (dotted). The 
adiabatic parameter is ^ = 0.043. 



equation, one would find that the system relaxes into the 
instantaneous ground state of the Hamiltonian. Because 
the two states have slightly different populations of the 
state |1) (the solid blue line does not approach the gray 
dashed line in the limit t — > 00), the importance of tak- 
ing into account the time dependence of the Hamiltonian 
can be tested experimentally by measuring cr^ after a 
sufficiently long time. This effect scales linearly with the 
adiabatic parameter (see section |ll| and is particularly 
pronounced if one is in the slightly non-adiabatic regime, 
which is the case in Fig. [T] 



B. Cooper— pair pumping 

In a Cooper pair sluice pTI - f25] [see Fig. [2] (a) for the 

design], the pumping of Cooper pairs is achieved by a 
periodic and adiabatic modulation of three system pa- 
rameters. The effective Josephson couplings Jl,r. can 
be controlled by magnetic fluxes through the left and 
right SQUID, respectively, while the charging energy 



Ecin - UgY 
gate voltage. 



is controlled via n„ 



Here, E, 



c 



2eVC, 



which is tuned by a 
where C is the total 
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FIG. 2: (Color online) (a) Schematic circuit of a Cooper pair 
sluice with a superconducting island (green) coupled to a su- 
perconducting circuit by two SQUIDs. The electrostatic po- 
tential of the island is controlled via a gate voltage, (b) The 
pulse sequences for the tunnel rates and the electrostatic po- 
tential. Throughout the paper, we use Ec/{2nh) — 21 GHz. 



capacitance of the island, and h is the number operator 
of the excess Cooper pairs on the superconducting island. 
Furthermore, the pumped charge depends on the phase 
difference (p of the two leads, which is fixed by embed- 
ding the design into a superconducting loop threaded by 
a fixed magnetic flux. In the following we only consider 
if = n/2. 

If the Josephson couplings are small compared to the 
charging energy, then a two level approximation is valid, 
where only the two eigenstates of the number operator 
with the lowest charging energies are considered. The 
Hamiltonian is then |22l 



H 

By 



-^{Bx(Jz 
Jl + Jr. cos ip, 

Ec{l-2ng). 



By ay 



(24) 



The pumping pulses are shown in Fig. [2] (b). Essen- 
tially, Hg is raised while the left Josephson junction is 
open, therefore the island gains one Cooper pair from 
the left lead. When Ug is lowered, the right junction is 
open and a Cooper pair is pushed into the right lead. We 
use a sufficiently slow pumping rate to assure good adi- 
abaticity [24| and without coupling to an environment, 
the system stays in the Floquet ground state if it was ini- 
tially prepared in this state. The resulting current ideally 
is one Cooper pair per pumping cycle. 

In an experiment, the Josephson junctions can not 
be closed exactly, and we assume a realistic ratio of 
J mini J max = 0.03. Therefore, in addition to the pumped 
current which depends on the pumping frequency, there 
will be a supercurrent which is independent of the pump- 
ing frequency [3S]. This different frequency dependence 
can be used to extract the pumped charge. 

In this system the main source of decoherence are 
fluctuations of the gate charge Ug, which according to 




FIG. 3: Panel (a) shows the energy splitting in the two-level 
approximation in units of 2-nh. Panel (b) shows \ {4>+ \ Oz \ 
(solid) and \\ (0+| |(^+) — {<\>-\ | (dashed), which re- 

late to the dissipation and dephasing rate, respectively. 



Eq. (24) results in a coupling to the environment de- 
scribed by the operator Eco^- Furthermore, we as- 
sume that the antisymmetrized environmental correla- 
tion spectrum \ [7(0;) — 7(— w)] = 70W increases linearly 
with frequency, as is typically the case for voltage fluctu- 
ations. We also assume an environment at thermal equi- 
librium 7(— w) — e~''"/^'^^'^^7(w), which leads to 7(w) = 



70^ 



1 — exp 



\kBT J 



as well as 7(0) = ^oksT/h, al- 
though dephasing is not important in a Cooper pair 
pump because the system starts in the ground state. 

It was already shown in [10| that a zero tempera- 
ture environment has no effect on ground state pump- 
ing. Here, we consider a flnite temperature environment, 
which is able to excite the system. In Fig. [s] (a), the 
level splitting of the Cooper pair sluice is plotted as a 
function of time for one pumping cycle. It is important 
to note that at times when the level splitting is small, 
i.e. when the system can be excited by a relatively low 
temperature environment, then the dissipative coupling 
to the gate voltage fluctuations has a pronounced max- 
imum, as is seen in panel (b) of Fig. |4] Therefore, we 
can expect that voltage fluctuations with temperature 
T = 200 mK « 4 GHz x2'Kh/kB already have a signifl- 
cant effect on the state of the system. 

Transient phenomena. The first few pumping cy- 
cles are not of importance in Cooper-pair pumping, as 
the experimentally measurable quantity is the averaged 
current after many cycles. However, in many other appli- 
cations such as adiabatic quantum gates, one is interested 
mainly in the behavior of the system within one period. 
Therefore, we discuss here some interesting effects which 
happen before the system approached the quasi station- 
ary state. Furthermore, this discussion might also be 
helpful in understanding the quasi stationary behavior 
studied below. 

In panel (a) of Fig. |4] we consider the system initially 
in the Floquet ground state and we plot the population 
of the Floquet ground state as a function of time for a 
fairly strong coupling to the dissipative bath for differ- 
ent temperatures. As predicted in |10) . the system does 
not get excited if the environment has zero temperature 
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FIG. 4: Population of the Floquet ground state during one 
pumping cycle plotted as a function of time in pico seconds. 
Initially the system is prepared in the Floquet ground state. 
In panel (a) we compare different temperatures ksT /{2-Kh) = 
4 (solid), 2 (dashed), (dotted) GHz at coupling to dissi- 
pative bathe (voltage fluctuations) of strength Eq^q — 0.1, 
while in panel (b) we compare different coupling strengths 
-Ec7o = 0.2 (solid), 0.05 (dashed), 0.02 (dotted) at tempera- 
ture kBT/{2Tvh) = 4 GHz. 



(dotted line). As expected, the higher the temperature, 
the more populated the Floquet excited state becomes, 
in particular at times when the level splitting is small. 

In panel (b) of Fig. |4] we plot the population of the 
Floquet ground state as a function of time for a temper- 
ature of about 200 mK and various environmental cou- 
pling strengths. If the coupling is strong (solid line), the 
system follows the instantaneous thermal mixture of the 
Floquet ground and excited states. For weaker couplings 
(dashed and dotted lines), the system gets less excited 
during the short time where the level splitting is small. 
Furthermore, the short time interval with dissipative cou- 
pling [see Fig|3](b)] is not sufficient for the system to re- 
turn to the ground state once the level splitting is large. 

Finally, in Fig. [5] we use even weaker coupling to the 
environment such that the system needs several pump- 
ing cycles to approach a quasi stationary state. As men- 
tioned in section [nij this quasi stationary state still has 
time dependent populations of the two Floquet states. 

Quasi stationary state. Now we assume that the 
system has already reached its quasi stationary state. In 
panel (a) of Fig. [6] we show the time dependence of the 
population of the Floquet ground state for various tem- 
peratures. Again, at times of a small level splitting the 
population of the Floquet excited state is increased, while 
at times of a large level splitting the system returns to- 
wards the Floquet ground state. 
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FIG. 5: Population of the Floquet ground state for the first 
(dotted) and second (dashed) pumping cycles, as well as for 
the quasi stationary state (solid). We use the parameters 
kBT/{2Tvh) = 4 GHz and £^70 = 0.01. 
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FIG. 6: Population of the Floquet ground state in the quasi 
stationary state. In panel (a) we compare different temper- 
atures fcsT/(27rfi) = 2 (dotted), 4 (dashed), 6 (soUd) GHz 
for a coupling to dissipative bath (voltage fluctuations) with 
Ec-yo = 0.1. In panel (b) we compare different coupling 
strengths E^^fo = 0.2 (solid), 0.02 (dashed), 0.002 (dotted) 
at temperature ksT / {2Tih) = 4 GHz. 



In panel (b), we compare the cases of different cou- 
pling strengths to the dissipative bath. If the coupling is 
strong (solid) the system closely follows a thermal mix- 
ture of the two Floquet states according to the instan- 
taneous level splitting. On the contrary, if coupling to 
the environment is weak (dotted) , it barely influence the 
system's state within one pumping cycle, which leads to 
a fairly constant population of each Floquet state. This 
is the limit which is described in Ref. [12j. The dashed 
line represents the intermediate situation. 

Current and pumped charge. The operator of the 
current through the right Josephson junction is given by 
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FIG. 7: (Color online) The expectation values of the observ- 
ables (Jx (black, solid), ay (red, dashed), and (blue, dot- 
ted) are plotted as a function of time. For panel (a) the 
density operator is calculated without using the secular ap- 
proximation, while for panel (b) the secular approximation 
was applied. The similarity of the two panels verifies the ac- 
curacy of the secular approximation when being interested 
in the state of the system. The small peak of at around 
J = 30 ps is responsible for the current and is shown in detail 
in panels (c) and (d). We see that this peak is altered by 
the secular approximation, which explains why this approxi- 
mation should be avoided when calculating the current. We 
use EqJo = 0.2 and fesT/(27rfi) = 4 GHz. For clarity, we 
use Jmin = to avoid supercurrent and to expect one trans- 
ferred Cooper pair per cycle (see main text). Indeed, panel (a) 
predicts 1.00 Cooper pairs, while performing the secular ap- 
proximation [panel (b)] leads to the incorrect value of 0.62 
transferred Cooper pairs. 



(see, e.g., Ref. [HI [52]) 

^ _ 2edH 

h dip 



ieJr. 



-lip 



(25) 



For = 7r/2 we obtain / = —^^a,j,. The expectation 
value of the current is rather small if the state of the 
system is close to the ground state and therefore small 
errors in the density matrix can lead to significant rel- 
ative errors in the current. In particular, it was argued 
in [ini [H] and discussed in detail in [55] that the secu- 
lar approximation leads to incorrect values for the cur- 
rent and therefore to the non-conservation of charge. In 
Fig. [7| we analyze the validity of the secular approxima- 
tion. We observe that the secular approximation is very 
good for the expectation values of all three observables, 
Uxi <7y, and fJz- Thus the quantum state of the system 
is well approximated. However, the small deviations in 
{ax) shown in Fig.[7]are sufficient to significantly modify 
the result for the pumped charge, when integrated over 
long time. 

We discuss the situation shown in Fig.[7|in more detail. 
To facilitate the qualitative understanding we assume, 
first, a vanishing residual Josephson coupling J,„in = 
and a coupling to the environment which is sufhciently 
strong for the state to closely follow the instantaneous 
thermal equilibrium. The temperature is of the order of 



the minimal energy gap, but much lower than the en- 
ergy gap at the beginning and the end of each half-cycle 
[see panel (a) of Fig. [S]. At the beginning of the cy- 
cle the system will therefore be in the ground state with 
no excess Cooper pair on the island. When the system 
moves through the first avoided level crossing, the en- 
ergy gap is small enough for the excited state to be pop- 
ulated. Nevertheless, after the first half-cycle is com- 
pleted, the energy splitting is large again and the system 
will again be in the ground state, which now has one 
excess Cooper pair on the island. Because Jji — Q dur- 
ing the the first half-cycle, this Cooper pair came from 
the left lead. With the same reasoning, the Cooper pair 
leaves the island to the right lead during the second half 
cycle. Therefore, in this parameter regime we should ob- 
tain exactly one transferred Cooper pair per pumping 
cycle, much like in the ground state pumping without 
coupling to an environment. Indeed, avoiding the secu- 
lar approximation we obtain exactly one Cooper pair per 
cycle [panel (a) of Fig. [t], while performing the secular 
approximation [panel (b)] leads to the incorrect value of 
0.62 transferred Cooper pairs. 

This failure of the secular approximation discussed 
above is not necessarily connected to the time depen- 
dence of the Hamiltonian. To illustrate this, consider a 
time independent Hamiltonian with = and Jr ^ 
and assume the system to be relaxing from some initial 
non-equilibrium state to the thermal equilibrium state. 
In this process the populations of the excited state and 
the ground state change in time. Since these states have 
different expectation values of charge on the island, some 
current must flow. However, for = only off-diagonal 
elements of the density operator in the energy basis are 
responsible for current and in the secular approximation 
the off-diagonal elements are never generated if they were 
not present in the initial state. This explains why the sec- 
ular approximation must fail to take into account contri- 
butions to the current which are induced by dissipative 
effects. 

To guarantee the correctness of the transferred charge, 
in the following we use the master equation Eq. (17) 



which was derived without the use of the secular approx- 
imation. 

In absence of dissipation the pumped charge Qp and 
the charge transferred due to the super current were 
found [231 US] to be given by 



Qp = -2e 1 - 2 



t/rr 



COS if 



Qs OC T/minSin(/3, 



(26) 
(27) 



where Imin = 2eJ,„in//i is the residual critical current 
of the system. In the strictly adiabatic limit, these for- 
mulas need not to be modified if the environment has 
zero temperature [TD]. Our choice of (^9 = 7r/2 then leads 
to exactly one pumped Cooper pair per cycle, and any 
deviation from that ideal value is because the finite tem- 
perature environment is able to excite the Cooper pair 
sluice. 
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FIG. 8: (Color online) The total current is plotted as a 
function of time. Panel (a) compares different tempera- 
tures kBT/(2-Kh) = (black, solid), 2 (red, dashed), 4 
(blue, dotted), 6 (gray, solid) GHz at fixed coupling strength 
5^70 = 0.05. Panel (b) compares different coupling strength 
Ec'yo — (black), 0.005 (red), 0.1 (blue) at fixed temperature 
ksT / (2nti) = 4 GHz. The insets highlight the second current 
peak. 

In Fig. [S] we plot the current in units of Cooper pairs 
per picosecond as a function of time, assuming that the 
system has already undergone enough cycles to be in the 
quasi-stationary state. In panel (a) we compare differ- 
ent temperatures and we find that as the temperature 
increases, the current through the sluice decreases. The 
reason is, that the Floquet-excited state carries supercur- 
rent in the opposite direction to that of the ground state. 
In panel (b), we compare different coupling strengths to 
the dissipative bath (voltage fluctuations), and similarly 
find that the peak current decreases with the coupling 
strength. However, the peak width increases with cou- 
pling strength such that the total transferred charge per 
cycle increases with the coupling strength (see panel (a) 
of Fig. |9]). This behavior can be explained qualitatively 
by looking at Fig. |6] (b). At the time of maximal current, 
the population of the ground state decreases with grow- 
ing coupling strength, therefore leading to a decreased 
peak current. But at other times, the population of the 
ground state increases with growing coupling strength, 
as does therefore the current, resulting in wider current 
peaks. 

Panel (a) of Fig. [9] shows the total transferred charge 
per pumping cycle Qtot = Qs+Qp plotted as a function of 
the coupling strength to the dissipative bath for different 
temperatures. Again, the charge decreases with tempera- 
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FIG. 9: (Color online) Panel (a) shows the total trans- 
ferred charge per pumping cycle in units of Cooper pairs 
as a function of the coupling strength EqJo to the envi- 
ronment for the temperatures ksT / {2Tih) =Q GHz, 2 GHz, 
4 GHz, 6 GHz (from top to bottom). The residual coupling 
is Jmin = O.OSJmax. Pauel (b) shows the pumped charge per 
cycle (solid curves), which is calculated from the difference of 
the total charges obtained by traversing the pumping cycle 
in opposite directions. For comparison (blue, dashed curves) 
we also show the transferred charge per cycle without resid- 
ual current Jmin = 0, which shows that the pumped charge is 
barely influenced by the residual coupling Jmin. 

ture, but it slightly increases with i?^7o. Panel (b) shows 
the corresponding pumped charge Qp, which is obtained 
from Qp = (Qtot — Qtot)/^: where Q^ot the transferred 
charge if the pumping cycle is traversed in the opposite 
direction which changes the sign of the pumped charge 
but leaves the supercurrent unchanged. As expected, the 
pumped charge decreases with temperature. The depen- 
dence on the coupling strength to the environment shows 
an increase with the growing coupling. This supports the 
conclusion of |10j that a strong coupling to the environ- 
ment might be helpful for Cooper pair pumping, and can 
be explained by the fact that at times when the level 
splitting is large, the system relaxes more effectively to 
the ground state. As elaborated before for Jmin = 0, for 
one pumped Cooper pair per cycle, it is sufficient if the 
system is in the ground state at the beginning of each 
half cycle. We further plotted the transferred charge for 
Jmin = which shows that the pumped charge is barely 
influenced by Jmin- For kBT/{2iTh) =6 GHz and strong 
coupling the pumped charge does not converge to one 
Cooper pair, but approaches a somewhat lower value. 
That is because a higher temperature environment is able 
to excite the system even at the end of each half cycle 
when the energy splitting is large. 



V. CONCLUSIONS 

In conclusion, we have explicitly taken into account the 
time dependence of an adiabatically and cyclically evolv- 
ing Hamiltonian to develope a master equation which is 
valid if the environment is Markovian. We found that 
decoherence takes place in the time dependent Floquet 
basis and that the decoherence rates depend on time 
through the instantaneous level splittings as well as the 
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instantaneous coupling to the noise. 

We first applied our theory to a spin in an adiabatically 
changing magnetic field, where we found that the spin 
does not relax into the instantaneous ground state, but 
the slightly difi^erent Floquet ground state which is very 
close to the superadiabatic ground state. We then tested 
our theory on the example of the Cooper pair sluice, 
where the time dependence of the decoherence rates is 
very pronounced. This results in new features at finite 
temperatures and modifies the pumped charge which is 
measurable in current experiments [HI US] . We further 
explain why the current has to be calculated without us- 



ing the secular approximation despite the fact that this 
approximation is very well justified if one is interested 
only in the state (density matrix) of the system. 
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